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ABSTRACT 

We apply the Wiener Hermite (WH) expansion to the non-linear evolution of the 
large-scale structure, and obtain an approximate expression for the matter power spec- 
trum in the full order of the expansion. This method allows us to expand any random 
function in terms of an orthonormal basis in the space of random functions in such a way 
that the first order of the expansion expresses Gaussian distribution, and others are the 
deviation from the Gaussianity. It is proved that the WH expansion is mathematically 
equivalent to the T-expansion approach in the renormalized perturbation theory (RPT). 
While exponential behavior in the high-fc limit has been proved for the mass density 
and velocity fluctuations of dark matter in the RPT, we prove the behavior again in 
the context of the WH expansion using the result of the standard perturbation theory 
(SPT). We propose a new approximate expression for the matter power spectrum which 
interpolates the low-k expression corresponding to the 1-loop level in SPT and the high- 
k expression obtained by taking a high-A; limit of the WH expansion. The validity of our 
prescription is specifically verified by comparing with the 2-loop solutions of SPT. The 
proposed power spectrum agrees with the result of iV-body simulation with accuracy 
better than 1% or 2% in a range of the baryon acoustic oscillation scales, where the wave 
number is about k = 0.2-0.4 /iMpc -1 at z = 0.5-3.0. This accuracy is comparable to or 
slightly less than the ones in the closure theory, the fractional difference of which from 
the iV-body result is within 1%. One merit of our method is that the computational 
time is very short because only single and double integrals are involved in our solution. 



1. Introduction 

Precise measurements of matter power spectrum in the large-scale structure are a powerful tool 
not only to investigate the details of the structure formation, but also to estimate the cosmological 
parameters. For example, the precise measurements of the baryon acoustic oscillation (BAO) in the 
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matter power spectrum observed by the Sloan Digital Sky Survey has emerged as a powerful tool 
to estimate cosmological parameters 

aaaaaaa 

la). Also the observation of cosmic shear 
in the near future is expected to give a useful constraint on the nature of dark energy. Obviously, 
proper understanding of the observed power spectrum becomes possible only if an accurate theo- 
retical prediction is available which requires a good understanding of the non-linear evolution of 
dark matter perturbation, the relation between dark matter and baryonic matter (bias effect) and 
the redshift distortion effect. There have been various studies and much progress on the theoretical 
calculations of the power spectrum, but it is still useful and required to have more accurate theoret- 
ical treatment. In this paper, we give a new approach to describe the non-linear evolution of dark 
matter. It is called the "Wiener Hermite (WH) expansion method" , where the stochastic nature of 
the cosmological density perturbation is manifestly used and the stochastic variables are expanded 
in terms of an orthonormal basis in the space of stochastic functions. The method was developed 
in the 1970s for application to turbulent theory in fluid dynamics and applied to cosmological tur- 
bulent theory by one of the authors of this paper. The method gives us a coupled equation at each 
perturbative order, even at the first order in such a way that the lower order quantities are modified 
by higher order quantities. This is totally different from the usual perturbation theory where lower 
order quantities are never influenced by higher order quantities. Thus, it gives us a prescription 
for the renormalization of higher order effects, and the precise meaning is described below. Each 
expansion coefficient has a clear statistical meaning; namely, the coefficients of the first, second 
and third terms in the expansion express the amplitude of Gaussianity, the skewness and kurtosis, 
respectively. Thus each term corresponds directly to an appropriate n-point correlation function. 

We mention here some details on the previous approaches relation to ours. It has been known 
for some time that the standard perturbation theory (SPT) of cosmological perturbation can be 
analytically solved in the Einstein-de Sitter universe in integral forms (ja; 10; II ; 12; H; 0; 0)- 
When it is considered up to the third order in SPT (1-loop level), the analytical predictions describe 
the nonlinearity well at sufficiently high redshifts (0; However, the predictions are still 

insufficient at the observable low redshifts (z = 0-3), and we need to consider further non-linear 
effects. Furthermore, it is computationally expensive to deal with the higher order corrections in 
SPT. Therefore, various modification of SPT have been proposed in the past. One of the main 
approach is the "Renormalized Perturbation Theory" (RPT; (|18l;ll9l;l20l)). where the basic equations 
for fluid describing matter perturbation are rewritten in a convenient compact form in order to use 
a diagrammatic technique developed in quantum field theory (12 ll ). Further modification have 
been considered, such as, e.g., the "Closure Theory" (1221: l23l). the "Time Renormalization Group" 



approach (J24j), and the 'T-expansion approach" using Multi-Point Propagators (|2a ; 
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28|). 



On the other hand, there is also an 



Many other new methods have also been studied (|29j; l30l : 
approach to the large-scale structure in the framework of the Lagrangian picture, called "Lagrangian 
Resummation Theory" (LRT; (0; 0; 0; 0)). 



It will be shown that our approach is mathematically equivalent to the T-expansion approach, 
but it still has the features described above and gives us a very convenient expression for the matter 
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power spectrum described below. In almost all modified perturbation theories, the resummation of 
nonlinear effects, which means the partial summation of the infinite order in SPT, is considered. 
This implies that any modified perturbative expansion methods should be described in the context 
of SPT. In this paper, we use only SPT, and prove again various properties of cosmological per- 
turbations, e.g., their behavior in the small-scale limit (high-A; limit) proved in the context of RPT 



iialiafctfD. 



Since low- A; solutions can be safely computed using SPT, the derivation of more precise solutions 
of cosmological perturbations by interpolating between the 1-loop results and the high-A; behavior 



(|19l ; l28l ) has been attempted. However, some arbitrariness have remained for this prescription. To 
resolve this problem, we propose a unique interpolation between the low- A; solutions and the high-A; 
ones by assuming that the higher order solutions in perturbation theory are well approximated by 
the ones in the high-A; limit. Then, we precisely compute only up to 1-loop level corrections in 
SPT and replace the higher order corrections with the ones calculated in the high-A; limit. In this 
way we obtain an approximate full power spectrum, and the power spectrum shows a very good 
agreement with iV-body results up to rather high-A; (about < 0.2-0.4 /iMpc^ 1 ) within 1 % or 2 % 
accuracy. 

This paper is organized as follows. In Section 2, we first explain the stochastic properties which 
should be satisfied by the density and velocity perturbations of dark matter. In Section 3, we briefly 
review SPT, which will be used later. Then the WH expansion technique is explained in our context 
in Section 4. The relationship between SPT and the WH expansion method is established there, 
and we also show the mathematical equivalence between the WH expansion and the T-expansion. In 
Section 5, we prove again the high-A; limit behavior of the cosmological perturbations in the context 
of SPT and propose an approximate full power spectrum, where the lower order corrections are 
calculated only up to 1-loop levels in SPT and the higher order corrections are replaced with the 
high-A; solutions. In Section 6, we compare our result with some other analytic predictions and 
A^-body simulations. We compute the two-point correlation function in Section 7. We summarize 
our work and discuss future works in Section 8. 



2. Stochastic Nature of Cosmological Perturbations 

After decoupling, baryon and dark matter fluctuations are tightly coupled by the gravitational 
force, and the evolution can then be described by pressureless fluid equations (continuity equation 
and Euler equation) with the Poisson equation for the Newton gravity. Thus our basic equations 
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are as follows (|15l): 



t + v. H = o, 

_ + v Vv = -V0, 

V 2 + Ac 2 = AnGp, (1) 

where v, and <fi denote the mass density, velocity and gravitational potential, respectively, and 
A is the cosmological constant. 

When we transform the spatial coordinates as x — > ax and redefine the velocity as v = dx + u, 
where a is the scale factor and u is the peculiar velocity, we can express Eq. dTJ) as 

+ 3W(r)p(r, x) + V • [p(r, x)u(r, x)] = 0, (2) 
9U(T ' X) + W(t)u(t, x) + [u • V] u(r, x) = -W(r, x), (3) 



V72J, 

„2 



^ + Ac^ = 4vrGp, (4) 



a 

where the conformal time r is defined as adr = dt, and the conformal Hubble parameter H is defined 
as T~L = aH, where H is the Hubble parameter. We further defined the cosmological gravitational 
potential as $ = (f> + ±H'x 2 . 

In the standard cosmological perturbation theory, physical quantities are decomposed into the 
background part and the perturbative part. The background part of the mass density p is defined 
as 

p(t)^(p(x,t)) = (p(0,t)). (5) 

where (• • • ) denotes the ensemble average and we used the translation symmetry of the ensemble 
average. On the other hand, the peculiar velocity u has no background part because of rotation 
symmetry in the average sense. Therefore, the perturbative part of the mass density and the 
peculiar velocity has the property that their ensemble average are zero by definition: 

(Sp) = (u) = 0. (6) 
Averaging the above set of equations, we obtain the following background equations. 

^ + mp = o, (7) 

&R 4vrG 2 1 . 2 2 

^ = -— pa 2 + -AcV, (8) 
Integrating Eq. (jHJ), we find the usual Friedman equation, 

^ + c 2£ = 8^ a 2- + l Ac 2 a2) (g) 
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where the integral constant K, is interpreted as the spatial curvature. 

By subtracting the background equations from Eq. (J2j) ([3]), we find our basic equations in 
Fourier space as follows; 

77^3 J ^3(2vr) 3 fe(k 1 + k 2 -k)a(k 1 ,k 2 )0(r,k 1 )<5(r,k 2 ), (10) 
9'(r,k) + H6(T,k) + -n m U 2 5(r,k) = - j ^3 J ^3(27r) 3 fe(k 1 +k 2 -k)/3(k 1 ,k 2 MT,k 1 )0(T,k 2 ) 

(11) 

where <5 = 5p/p and = <9jU* denotes the divergence of velocity, and 6d denotes the three- 
dimensional Dirac delta distribution. We neglected the vorticity w = V x u because the vorticity is 
zero if its initial value is zero, and even if its initial value is non-zero, it decays due to the expansion 
of the universe. The functions 

_ (ki + k 2 ) ■ ki 
a(ki,k 2 ) = ^ , (12) 

an , x _ |kj +k 2 | 2 (ki -k 2 ) 

P^M) = ^I2l2 > (13) 

encode the nonlinearity of the evolution and satisfy the conditions 

a(k, -k) = /3(k, -k) = 0. (14) 



Note that this decomposition between background and perturbation is exact in Newtonian 
gravity, namely there are no backreaction terms generated from the ensemble average in Newto- 
nian gravity, and the perturbative parts 5 and 9 obeying Eqs. (fTUj) and (fTT|) naturally satisfy the 
stochastic condition in Eq. ([6]): (5) = {9) = 0. More specifically, when the solutions for nonlinear 
equations are expanded perturbatively, the property that their ensemble averages are zero is not 
guaranteed in general, and thus a redefinition of the perturbation variables such as 5 —> 5 — {5) is 
necessary. In the case of cosmological SPT, the average, such as (<5(k)) being proportional to <5/)(k), 
is interpreted as the vacuum bubble diagram in the diagrammatical picture and contributes only in 
infinitely large scales. This means we really need to redefine of the cosmological background. How- 
ever, we do not need to consider this prescription for perturbative variables in Eq. (flQl) and (fTT|) . 
This is a special feature of Newtonian gravity. It will be interesting to see how the backreaction 
look like in the case of general relativistic gravity in our approach. 



3. Review of Standard Perturbation Theory 

We here explain SPT very briefly which will be used later. The solutions in the case of an 
Einstein de Sitter universe, where O m = 1 and = 0, can be described by analytically integral 
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forms in SPT. More explicitly, the solution may be written in the following perturbative form, 

oo oo 

5(z,k) = £y<y n (k), 9(z,-k) = -Hj2a n 6 n (k), (15) 



n=l 



n=l 



where the scale factor a is a growing mode solution in the linearized theory. When the scale factor 
a is small, the series are dominated by their first term, that is, by linearized theory. The relation 
between the time-independent coefficients 5i(k) and #i(k) is shown from the continuity equation 
(fit))) as <5i(k) = #i(k) = <5i(k), and the time-independent linear power spectrum for <5z,(k) is defined 
as 

(^(k)^(k')) = (2vr) 3 5 D (k + k')P L (k), 



(16) 



where the amplitude of the wave vector is expressed as k = |k|. 

Then, the coefficients 5 n (k) and 9 n (k) are described as follows, 





d 3 q n 


(2vr) 3 


(2vr) 3 


d 3 Ql 


d 3 q n 


(2tt) 3 


(2tt) 3 



(2vr) 3 J jD (k - qin)i ? n(qi, • • • , qn)<5i(qi) • • • S L (q. n ) 
{2-jr) 3 5 D (k - qi„)G n (qi, . . . , q n )*i(qi) • • • S L {q n ) 



(17) 



(18) 



where qi n = qi + q2 + • • • + q n and F n and G n are completely symmetrized functions for the 
wave vectors {qi, q 2 , • • • , qni- The functions F n and G n are constructed according to the following 



recursion relations (n > 1) (jlOl : Il5l): 



F n+1 (q l5 . . . , q n+1 ) - ^ (2n - 5)n 



m=l 



(2n + 3)a(ki,k 2 )F n+ i_ m (q m+ i, . . . ,q„ + i) 

+ 2/3(ki,k2)G n +i_ m (q m +i, . . . , q n +i) 



(19) 



r, i \ (qi , • • • , qri 

G n+1 ( qi , . . . , q n+1 ) - ^ 

m=l 



3a(ki,k 2 )F rt+ i_ m (q m+ i, . . . ,q n +i) 

+ (2n + 2)/3(ki, k 2 )G n+ i_ m (q m+ i, 



i qn+i, 
(20) 



where k x = q x H h q m , k 2 = q m+ i H h q n+ i, and F x = d = 1. 

For n = 1, we have 

5 lki k 2 (k x k 2 \ 2(ki-k 2 ) 



F 2 ki,k 2 = - + - 



+ + = 



7 2 fcife \k 2 hj 7 k\k\ ' 



(21) 
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3 1 ki • k 2 fh k 2 \ 4 (ki • k 



\2 



^•*>-j + i^5 + £J + r=isr (22> 

The stochastic property (5) = (0) = is specifically shown from these solutions. When we 
consider the average of Eq. ([T7]). we only have to consider both coefficients 5 n and n due to the 
linearity of the ensemble average. The ensemble average for 5 n is 



<*n(k)> = J |^---^(2vr) 3 ( 5 D (k-q ln ) J F„(qi,...,qn)(5L(qi)---5L(qn)} 
d 3 qi d 3 q n 



(2vr) 3 fe(k - qi„)F n ( qi , . . . , q n )(27r) 2 5 D (q ln ) J B(q 1 , . . . ,q n ), (23) 



(2tt) 3 ' ' ' (2vr) s 

where -B(qi, . . . , q n ) is defined as 

(^(qi) • • • <5 L (q n )) = (2vr) 3 5 D (q ln ) J B(q 1 , . . . ,q n ). (24) 

When the function F n in Eq. (fl~9l) substitute into Eq. (|23|) . ki + k 2 = qi n = is satisfied due to 
the Dirac delta function in Eq. (|23p . and the functions a(ki,k2) and /3(ki,k2) become zero from 
Eq. (|14p . Then, it is shown that the function F n in Eq. ()23[) becomes zero, and {5) = 0. A similar 
analysis can be applied to 9. 

Note that this stochastic property is independent of the initial conditions of 8l, that is, the 
initial condition can have primordial non-Gaussianity. 

4. The Wiener Hermite Expansion 

Now, we explain our expansion method for 5 and 0. Our expansion scheme should satisfy the 
following two properties. First, it is known observationally and theoretically that the cosmological 
perturbations in the universe have a nearly Gaussian distribution. Thus the first order in our 
expansion should express the Gaussian distribution. Second, the expansion scheme should respect 
the stochastic condition of the cosmological perturbations, (5) = (9) = 0. Based on these two 
conditions, we adopt the WH expansion as our expansion method. 

4.1. Definition of the Wiener Hermite expansion 

In the WH expansion, the perturbation variables 5 and 9 are expanded as follows, 
d 3 pi f d 3 p r 



5(z,k) = pJ^.-. J ^(2vr) 3 5 D (k-p lr )4 ) H (z,p 1 ,...,p r )^ W (pi,..., Pr 



(25) 



' k ) = E/ |^T" / ^(2^ D (k- Plr )< ) H (z,p 1 ,..., Pr ) J ffW( Pl ,..., Pr ), (26) 
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where the functions {r = 1, 2, . . . } are the stochastic variables. The first-order is a white 
noise function which satisfies the Gaussian distribution, 

( J ff( 1 )(k 1 ) J ff( 1 )(k 2 )) = (2^) 3 Mk 1 + k 2 ), (27) 
and we further define higher order bases 

ffW { r = 2,3,4,...} in the expansion as follows: 
H^{k u k 2 ) = ^ (1) (k 1 )F( 1 )(k 2 ) - (2vr) 3 5 D (k 1 +k 2 ), 

fl-( 3 )(k l5 k 2 ,k 3 ) =^ 1 )(k 1 ) J ff( 1 )(k 2 )F«(k 3 ) -F( 1 )(k 1 )(2vr) 3 fe(k 2 + k 3 ), 

- J ff( 1 )(k 2 )(2vr) 3 5 D (k 1 + k 3 ) - J ff( 1 )(k 3 )(2 7 r) 3 ( ) D (k 1 + k 2 ), 
fl-( 4 )(k l5 k 2 ,k 3 ,k 4 ) = ... . (28) 

Thus they become an orthonormal basis in the space of stochastic functions, where the ensemble 
average of is clearly zero. 

( J ffW(k a ...k r ))=0 J {r = l, 2, 3, ...}. (29) 

(i2-W(k 1 ...k r )iT«(k / 1 ...k / s )} = 0, {r^s}. (30) 

Therefore, the stochastic property in Eq. (|6|) is satisfied by the definition of the WH expansion. 

The coefficients of the WH expansion are derived by averaging 5 and 6 after multiplying the 
stochastic variable : 

(S(k)^(-'k 1 , -K)) = (2vr) 3 5 D (k - ki r )r!<S# H (.z,ki, . . . ,k r ), 
^(k^t-kx,...,-^)) = (2irf5 D (k -k^rW^zM, ■ ■ ■ ,K)- (31) 



The power spectrum is described in the WH expansion by 

'/''/' I f d 3 p r rr (r+l)/_. ,„ _ „ _ ,,2 



P{z,k) = ]T(r + l)! J "^•••/ ^^[*w£ 1) («. k -Pir,Pi 5 --- J Pr 

oo 

= E P WH 1) (-^)> (32) 

r=0 



where the contribution of the (r + l)th order in the WH expansion to the power spectrum P A 



(r+l) 
WH 



is defined as 



P&%,k) = {r + l)l J | ^[«jSS 1) (*,k-p lr) p 1 ,... ) p r )] 2 . (33) 
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4.2. Relation between the Standard PT and Wiener Hermite Expansion 

Since the solutions for any order of SPT are given in an Einstein de Sitter universe analytically, 
any new expansion for 5 and 9 must be described in the context of SPT. For the linear order, we 
assume 

6 L (k) = 4 1} (fc)ff«(k), P L (k) = [6?(k)\ 2 , (34) 

where the superscript and script indices for 5^ denote the order of the WH expansion and SPT, 
respectively. It is straightforward to include the intrinsic non-Gaussianities as higher order contri- 
butions in the WH expansion. From now on, we express 5^ (k) — > Substituting the solutions 
in the SPT Eqs. (fT5|) . (fTTl) . and (fl8|) into Eq. (J3T1) . we find the general relation of the solutions 
between SPT and the WH expansion as follows: 

oo 

x( r+1 )f„ lr. 1^ ,\ — V s „2n+r+l jr(H-l) 



(z,k!,...,k r+1 ) = ^a 2 -+ r+l 4;+ 1 r +i( k i---x+i), 

n=0 

OO 

{zM,-.-,K+i) = -Hj2a 2n+r+l 4n+%iQ<U---*r+i), (35) 



n=0 



where 



fclUki, • • • , kr+i) = (r + 1)! (2 " 2 n+|)? )! (2 " + l)ll6L{h) ■ ■ ■ 6L{K+l) 

d 3 pi d 3 p n 



/d pi d p n 
(2vr)3 ' ' ' (2^ F2n+r+1 ^ kl ' ' ' ' ' kr+1 ' Pl ' ~ Pl '- • • ' Pn ' -Pn) P L(Pl) ■ ■ ■ PhiPn) 

4^ +1 (ki, . . . , k r+1 ) ^ (r ^ 1) , (2 " 2 ^t)! )! (2n + 1)ll5L{kl) ■ ■ ■ 5Ukr+l) 

/d^ i p~^ d^ f p 
(27^ ' ' ' T2n^ G2n+r+1 Q SLl > ' ' ' ' kr+1 ' Pl > ~ Pl ' • • • ' Pn ' -Vk)Pl{pi) • • • Phipn) 

(36) 

This expression means that the density fluctuation with order (r + 1) in the WH expansion is the 
sum of all the density fluctuations with order (2n + r + 1) in SPT. 

We can understand the relation between SPT and the WH expansion through a diagrammatical 
representation, where there is no non-dimensional coupling constant and the order of the loop is 
determined by the order of the linear power spectrum Pi(k) } that is, the n-loop contributions 
contains the terms proportional to (Pl) u+1 - Each order of the WH expansion includes all vertex 
loop contributions which come from the 5^ coefficients themselves. The order of the vertex loop 
is expressed by n (n > 0). On the other hand, the loop contributions from irreducible diagrams, 
where the loop order is expressed by r (r > 0), arise only after calculating the power spectrum in 



Eqs. J32D and $2fr (see d23 123) for details). 



-10- 



4.3. Relation between the T-expansion and Wiener Hermite expansion 



The relation given in Eq. (|3ip corresponds to Eq. (17) in (|25l ). That is, the WH expansion 
method coincides with the T-expansion approach: 

r( r+1 \z, ki, . . . , k r+1 ) = (z, k x , . . . , k^o/Ofcfa) • • • feCfer+O). (37) 

For r = in Eqs. ([35]) ([36]) and ([37]) . we have 
r^(k) = 5^ H (z,k)/5 L (k) 

d 3 pi d 3 p r 



5> 2 ™ +1 (2n + l)!! J ^EL . . . ±^p 2n+l{K pi; _ pi , . . . , p „, - Pn )P L ( Pl ) ■ ■ ■ P L ( Pn ). 

(38) 



n=0 



Thus, the first order of the WH expansion, <^ H , corresponds to the propagator in RPT. Further- 
more, ^r+i 1 ^ denotes the irreducible diagrams, and this is expressed in the T-expansion as 

r&^C*, ki, . . . , K+i) = a^S^Oct, k r+1 )/ (5 L (h) . . . 5 L (k r+1 )) . (39) 
Note that since we focus only on SPT, we consider only growing solutions, unlike RPT. 



5. Behavior of the solutions in the high-A; limit 

Although the WH expansion is defined and the interpretation is physically and mathematically 
useful for understanding the nonlinear evolution of dark matter, we need to truncate the expansion 
at some order, most probably at a lower order such as r = 2 or 3 in order to perform the actual 
calculation. However, the validity of the truncation is not guaranteed immediately. Furthermore, 
the computational difficulties for the power spectrum increase very rapidly when we increase the 
order of truncation. In order to resolve these difficulties, we propose in this section an approximate 
semi-analytic expression for the full power spectrum including all order in the WH expansion by 
adopting the following assumption: The high-order solutions in SPT become dominant in the high-fc 
limit. Therefore, they are approximated well enough by the ones in the high-/c limit. Here, we show 



the exponential behavior of the solutions using SPT, which have been proved in RPT (|la : 1 13 : 1271 ) 



5.1. Functions F n and G n in the High-fc Limit 

We prove that the functions F r+n+ i(k, ki, . . . , k r , pi, . . . , p n ) and G r+n+ i(k, ki, . . . , k r , pi, . . . , p n ) 
take the following expression in the high-A; limit: 

(r + 1)! 

F r+n+ i(k,ki, . . . ,k r ,pi, . . . ,p n ) -> - — ■ — -F r+ i(k,ki, . . . ,k r )7(p x ) . . -7(p n ), 

(r + n + 1)! 

(r + l) 1 

+n+i (k, ki , . . . , k r , pi , . . . , p n ) — y - — ■ ■ — — rrG r +i(k, ki , . . . , k r )'y(pi) . . . / y(p n ), (40) 

(r + n + lj! 
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with 



Here, we define the high- A; limit as 



7(P) 



p_k 

p2 



|k| > {|pi|, i = 1,2,... ,n}. 



(41) 



(42) 



From now on, we shall simplify the notations: F r+n+ i(k, ki, . . . , k r , pi, . . . , p n ) = F r+n+ i(k, k r , p n ) 
and G r+n+ i(k, ki, . . . , k r , pi, . . . , p n ) = G r+ra +i(k, k r ,p n ). 

We prove this by induction in n as follows. For n = 0, Eq. (|4U|) is clearly satisfied. For some 
n, we assume 

(r + 1)! 



-^r+n(k, k r , Pn— l) 
Cr+n(k, k r , Pn— l) 



(r + n)\ 
(r + 1)! 



(r + n)\ 

Then we show that the n + 1 order satisfies the same limit. 



F r+ i(k,k r )7(pi) . ..7(p n _i), 
G r+ i(k, k r ) 7 (pi) . . . 7(p n -i). 



(43) 



The functions F r+n+ i(k, k r , p„) and G r + n +i(k, k r , p n ) are given by Eqs. (fT9|) and (i20l) . Then, 
let us examine which terms become dominant in the high-fc limit in these recursion relations. 
From Eq. (|42|) . we keep only terms with scale dependence as (k/pi) ■ ■ ■ (k/p n ). Then, we have the 
terms proportional to F r+ „(k, k i; p n ), G r+n (k, k t , p„) for i < r, and j(p n )F r+n (k, k r , p n _i) and 
7(pn)G r + n (k, k r , Pn-i) in the high-/c limit. This means that the recursion relation for F r+n+ \ (k, k r , p r 
in the high-/c limit becomes: 



-^r+n+l(k, k r , Pn 

(2(r + n) + 3) 



1 



(2(r + n) + 5)(r + n) 
C(r, m) 



+ 4 



C(n + r + 1, m) 

r 

+ E G 

m=l 

C(r, m) 



C(n + r + 1, m 
+ (2(r + n) +3) 



+ 4 



^ ] G m (k m )a(ki m ,k + k 

(m+l)r + Pln)-Pr+n+l— m(k, k m _|_i, . . . , k r , Pn) 
r+n+1— m(k, k m +l, . . . , k r , p n )fi(k + k( m _|_;Q r + Pin; ki m ) F m (]si rr , 

G m (k m )/3(k lm ,k + k 

(m+l)r "I - Pln)C r _|_n+l— m(k, k m +l , • ■ ■ , k r , p n ) 

a(Pn, k + ki r + p 1 („_ 1) )F r+n (k, k r , p n _i) 

(44) 



m=l 

n 



r + n+1 



— j /3(p n ,k + k lr + pi( n _i))G r+ „(k,k r ,p„_i)|, 



r + n + 

where we denote k( m+1 ) r = k m+ i + • • • + k r and Pi( ra -i) = Pi + • • • + Pn-i, and define as 



C(n, r) 



r!(n — r)! 



(45) 
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Furthermore, the behavior of a and /3 in the high-A: limit is 



a(ki m , k + k( m+1 ) r + p ln ) 
a(k + k( m+ i) r + Pin, ki m ) 
/3(ki m , k + k( m+1 ) r + pi„) 
a(p n ,k + ki r +Pi(„_i)) 

/3(p n ,k + ki r +pi( n _i)) 



a(ki m ,k + k( m+1 ) r ), 
a(k + k( m+ i) r ,ki m ), 
/3(ki m , k + k( m+1 ) r ), 

7(Pn), 

7(Pn) 



and Eq (I44[) becomes 

^r+ra+l(k, k r , 

x |(2(r + n) +3) 
+ 2 



(r + 1)! 

(2(r + n) + 5)(r + n) \{r + n+l)\ 



7(Pl) ■ ■ ■ 7(Pn) 



l-m(lm+l > • • • > Qr+l) 

m=l 
r 

J]G m (q m )/3(k 1 ,k 2 )G r+ 

m=l 

+ (2(r + n) + 3)nF r+ i (k, k r ) + 2nG r+1 (k, k r ) L 



where {qi,.. . ,q r +l} = {k, ki, . . . , k r } and ki = q lm , k 2 = q( m+ i)( r+ i). 
From Eq. (I19D and Eq. (I20p . we can show the following relations. 



^ G m (q m )a(ki, k2)F r+ i_ m (q m+ i, . . . , q r +i 

r?l=l 

G m (q m )/3(ki, k2)G r +i- m (q m +i, . . . ,q r +i) 



m=l 



(46) 



(r + l)F r+ i (k, k r ) - G r+ i (k, k r ) (47) 



-- (3F r+1 (k, k r ) - (2r + 3)G r+1 (k, k r )) 



(48) 



Substituting Eq. (I47p and Eq. (j48|) into Eq. (146p . we can finally derive the following relation 
in the high-/c limit, 



(r + 1)' 

F r+n+1 (k,k r ,p n ) -)• - — ■ — — 7(pi)...7(p n )F r+ i(k,k r ). 

(r + n + 1J! 

Similarly, for G r + n +i we can show the following relation 

(r + 1)! 

G r+n+ i(k,k r , Pn ) -> - — ■ — — 7(p 1 )...7(p„)G r+ i(k,k r ). 

(r + n + lj! 



(49) 



(50) 



This ends the proof. 
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5.2. Power Spectrum in the High-fc Limit 

We calculate the coefficients of the WH expansion in the high-A; limit, 



2n+r+l c(r+l) 



n=0 

oo 



J] a 2 "+ r+1 <5 L (|k - k lr |)^(fci) . . . 5 L (k r 

n=0 



x F r+ i(k - ki r ,k r ) 



2"n! 



"6vr 2 



dpPi(p) 



exp 
exp 



2 



-J 5 L (z, |k - ki r |)5 L (z, fei) . . . (5 L (z,/c r )F r+ i(k - k ir ,k r ) 
^^(^k-kx^kx,...,^), (51) 



where we have used Eqs. ([36]) . (f4T)j) . and (2n)!! = 2 n n!. We define cx 2 as 

a v = J dpP L (z,p). 



(52) 



Note that we define the z-dependent quantities such as 5i{z,k) and Pi(z,k) by multiplying the 
scale factor a, but we assume that the scale factor can be replaced by the growth factor D(z) 
in the general cosmological models, for which S1a ^ 0: 5l{z, k) = aSiik) — > D(z)SrXk) and 
Pl(z,p) = cl 2 Pl(p) — > D 2 Pl(p). This relation in Eq. (|5ip is equivalent to Eq. (42) in (25) from the 
relation between and r tre e in Eq. (|39|) . 



Then, we describe the full power spectrum in the high-/c limit as 



P(z,k)^e W (-k^J2 P lrr +1) ^k), 



(53) 



r=0 



where 

p£ +1 >(*,fc) = (r + l)! 

= (r + 1)! 



(2^)3 
d 3 fci 
(27) 



3 ' ' ' 



d 3 &r rr(r+l)/ iii i \i2 

7^33 K+i ; (2,k-ki r ,ki,...,k r )j 

d?k 

^ [F r+1 (k - k lr , k r )] 2 Pl(z, |k - k lr \)P L (z, h)... P L (z, k r ). 

(54) 



(2nf 



P 1TT includes the contributions from all the irreducible diagrams. 
Here, we take further high-/c limits in Eq. (|54p . 

k»{|ki|, % = l,...,r}. 



(55) 
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Using Eq. ([30)1. we show 

^^^-P L (z,k). (56) 

Note that the limit of k ^> k r (r > 1) is equal to the approximation that the most effective region 
of k r to each integral is around k r — > 0. However, since k r have integral range of < k r < oo, 
there necessarily exist the case of k r ~ k in the integral. For ki ~ k 3> k r (k r > 2), we have 

d?k\ f d 3 k r 



(57) 



where we define k as k — ki r ~ k — ki = q. This is equal to the case of the high- A; limit. The same 
analysis is applied to arbitrary k r . Therefore, the factor (r + 1) is multiplied in Eq. (|56|) . 

From Eq. ([56)) and Eq. ([53)) . we finally derive 

00 (k 2 ^Y 



P(z, k) P L (z, k) exp(-k 2 a 2 v ) 



r=0 



= P L (z,k). (58) 
Surprisingly, the solutions in the high-A; limit cancel each other, and the full power spectrum reduces 



to the linear power spectrum. This fact is well known in the 1-loop level of SPT (jl2j ). but it is 
interesting that this cancellation also applies for the dominant terms in the high-/c limit in the full 
power spectrum. Of course, it is really not that the full power spectrum becomes the linear power 
spectrum in the high-A; limit, because we have chosen only dominant terms in the high-/c limit in 
the proof of Eq. (|40p and Eq. ()56p . and the sub leading terms can also affect the full power spectrum 
even in the high-A; limit. This result implies that the nonlinear corrections for the power spectrum 
generally tend to cancel each other, and result in small corrections as specifically known for 1 and 
2-loop cases of SPT. 



5.3. Approximate Full Power Spectrum 

We now propose an appropriate interpolation between the low-A; and high-A; solutions. The 
low-A-solutions are the 1-loop solutions in SPT, while the high-A; solutions are given by Eq. (I53p 
derived in the previous subsection. 

In order to have an expression applicable for the case of r = 0, we use Eq. (j4Qj) in the following 

3' 

F 2n+ i(k, pi, -pi, . . . , p n , -p n ) -> — — — — F 3 (k, p, -p)7(p 2 )7(-P2) • • • 7(Pnh(-Pn)- (59) 

(2n + 1)! 
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Then, for n > 1, we show 
(i) /1N 2(2n + l)H 



J 2n+1 



(k) 



(2n + 1)! 

4 1} (*) 2 



2 n n\ 



6L{k)3 J |^^(k,Pi,-pi)PL(pi)] [| ^7(P)7(-P)^(P) 

A; 2 f "' ' 

-^2 J dP P Up) 



n-l 



where we have denoted the 1-loop correction term in SPT as 

d 3 p 



5 w {k) = 3 s L{k ) J jUL F3 (k,p,- P )p L (p). 



(60) 



(61) 



Then, we derive the approximate solution of <5^ H as 



n=0 

-> 5 L {z,k) + 5 ( i\z,k) 



-k 2 al 



S L {z,k) 



25^\z,k) 
Vol 



exp 



E 

n=l 



fcV 2 



1 / fcV. 



2^2 \ " 



(62) 



where we have used the general growth factor D instead of the scale factor a. The contribution to 
the power spectrum P^l in Eq. (f3"3"j) is 



. . 24 1) (^,A:) 



(63) 



For low-A;, we can expand the exponential term as e fc2<j2 / 2 ~ 1 — /c 2 <7 2 /2, and this lead the 1-loop 
correction, 

6^ u (z,k)^5 L (z,k) + 5i 1 \z,k). (64) 

While, for the high-A; limit, <5^ H becomes coincident with i^e - ^ "/ 2 due to the good convergence 
of 5 L + 25^/(^)^0. 

Next, for r > 1, we use the approximation of Eq. f)53[) . Here, we further approximate P-j^ 1 ^ 
because it is expensive to compute the terms in the case of (r > 2) due to their large multiple 
integrals. Using the following approximation from Eq. (|40p . 



F r+ i(k - ki r ,k] 



2! 



, . . . , a. r ) 



(r + 1) 



■ 7 (k 2 )...,7(kr)i ? 2(k-k 1 ,k 1 ), 



(65) 
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we derive the approximate solution of P^n^ as, 



P^^fcj^exp^ 2 ) (r + 1)! 



exp (-k 2 a 2 ) (r + 1)! 



d 3 ki d 3 k r 



S; r +i (z,k - ki r ,k] 



(2vr) 3 ' ' ' (2tt) 3 L 
2! 2! 



, . . . , jv r j 



(r + 1)! (r + 1)! 



r+1 



(2^)3 



ex P (-^)^fe^i(fc 2 



fc 2 cr^ r! 



where we multiply in the factor (r + l)/2 for the same reason as in Eq. 
another 1-loop correction term in the SPT as, 

d?p 



(66) 

We have denoted 



P 22 (z,k)=2 J AP_[F 2 (k-p,p)] 2 P L (z,\k-p\)P L (z,p). 



(67) 



i 2 2 

Indeed, if for r = 1 the exponential factor is expanded as e "* ~ 1, the 1-loop correction, P 22 , is 
reproduced in Eq. (f66|) . 

Finally, we achieve the approximate full power spectrum, 



P( Z ,k) = ^ P WH\^ k ) 
r=0 

24 1) (^,fc) 



Paf(+£0 



5 L (z, k) 



$L(z,k) 



k 2 a% 

2S^\z,k) 
k 2 ol 



exp 



exp 



k 2 al 



k 2 al 



+ exp i—k 2 



2 2 \ P22(z,k) ^ 1 2 2 



r=l 



+ ^^[l-exp(-^)]. (68) 



This is the main result of this paper. This gives an appropriate interpolation between the low-fc so- 
lutions and the high-A; limit ones. We can derive the approximate solutions of each order in the WH 
expansion and SPT using approximation such as Eq. (|59p and Eq. (|65p . and therefore we call this 
method as the "Approximate Full Wiener Hermite (AFWH)" expansion method or "Approximate 
Full Perturbation Theory" . We can easily compute the solution of Eq. (|68p numerically, because 
the solution has only single or double integrals. 



6. Comparison with Other Analytic Predictions and iV-body Simulations 

We compare the approximate full power spectrum, P ap , in Eq. (|68p with some other analytic 
predictions and iV-body simulations. We mainly use iV-body simulations presented in I22I ). but 
in Sec. 16.31 we also use other iV-body results with higher resolutions in (|36l). It is plotted for the 
cosmological models with the Wilkinson Microwave Anisotropy Probe (WMAP) five year ([37I ). 
cosmological parameters: Vl m = 0.279, ft A = 0.721, Sl h = 0.046, h = 0.701, n s = 0.96, a s = 0.817. 
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k[h/Mpc] 



Fig. 1. — Comparison between P^h in Eq. ([63]) (red line) and [PwiilReg m Eq. ([7U|) (black line) 



(i) 



for P^Ij/Pj, at z = 1. The fractional difference, (P^h ~~ [PwnlRegVPLj is also plotted. 



(i) 



The iV-body simulation data with low-resolutions and high-resolutions in (|22l ) and (|36l ) were 
created by a public iV-body code GADGET2 (|38l ). Their initial conditions were generated by the 
2LPT code flag) at 

^ini — 31 and z- m i — 99, respectively. While the A^-body simulations with 
low-resolutions were computed with a cubic box of size l/i _1 Gpc containing 512 3 particles, the 
TV-body results with high resolutions, called Lll-Nll and L12-N11, contain 2048 3 particles and 
were computed by combining the results with different box sizes, 2,048 /i _1 Mpc and 4096 ^ _1 Mpc. 
Although the realization of the simulations with high resolutions is only 1, the simulations with 
low resolutions have the output data of 30 independent realizations, and consider the correction 
of the finite-mode sampling by (J40j). Therefore, the size of each error bar for A^-body results with 
low-resolutions becomes hard to see visually. For details of the A^-body simulation used in this 
paper, see Taruya et al. (22) and Okamura et al. |35l). 



6.1. Comparison with Other Analytical Predictions: 1-loop Level 



(|28l ) proposed a simple scheme to interpolate between the \aw-k and high-A; solutions, based on 
the T-expansion method. In the scheme, the solutions are regularized so that the low-A; solutions 
become the ones in SPT and the high- A; solutions become Eq. f)53[) . We call this scheme "regularized 
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r-expansion" . From Eqs. (20) (51) in (|28l ? ), its 1-loop solution for the power spectrum is given 

by 

-Fkeg = exp(-k 2 a 2 ) 



k 2 n 2 \ 2 

5l + S 3 + ^5 l ) +P22 



(69) 



Since the WH expansion and the T-expansion are completely equivalent to each other, we can 
understand this solution as the truncation up to the second order in the WH expansion and write 
^WH us i n S the regularized T-expansion as 



(i) 



WH 



Reg 



exp 



-k 2 a 2 ) 



h + $3 + 



k 2 a 2 



n 2 



(70) 



The difference between Eq. (|63p and Eq. (|70p is manner of interpolating between the solutions. 
Since the regularized T-expansion method is a heuristic scheme, we have used the approximation 
of Eq. (|551> . 



When we ignore the contributions of Pmt (r > 2), we derive the solution of the approximate 



WH expansion corresponding to the regularized T-expansion in Eq. (|69p as 



pW , p(2) 



h(z,k) d , 9 _ 9 [ exp 



k 2 u 2 



k 2 al 



+ eM-k 2 (7 2 v )P 2 2(z,k). (71) 



We plot these two solutions, P$&/Pl and [PwlilRcg/ Pl, and their fractional difference, (-P^H - 
[-PvraJReg) / Pl, at z = 1 in Fig.HJ On BAO scales (~ 0.2/iMpc _1 ), the fractional difference is within 
1 %. At high k, our solution becomes slightly larger than the regularized T-expansion. However, 
there is no means of investigating which results are more accurate in detail, because on such scales 
the amplitude of Pwjj decays enough due to the exponential factor to not largely contribute to the 
full power spectrum. 



In Fig. [2J we plot the various analytic solutions with the 1-loop level corrections and A^-body 
simulation result (blue dashed: SPT, green dashed: Regularized T-expansion, black solid: LRT, 



orange solid: P^L + P^ii in Eq. ([71]) : red solid: AFWH in Eq. ([68]) ; and black symbols: A-body 



,(2) 



result) at z = 1. □ To easily see the BAO, we plot the ratio of power spectrum to a smooth reference 
spectrum, P(k)/P nw (k), where the function P nw (k) is the linear power spectrum calculated from the 
smoothed transfer function neglecting the BAO feature in (|4ll ). To investigate the agreement with 



1 The power spectra of SPT and LRT (|32f ) are given by 

Piioop = Pl +P13 +P22, (72) 

P Lag = exp (-k 2 (rl) {P L + P13 + P22 + k 2 (jlP L ) , (73) 

where we denote as P13 = 2<5_l<53 1 '. Although LRT is very similar to the regularized T-expansion and our result, the 
complete correspondence (e.g., the origin of the exponential factor) is not trivial. 
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Fig. 2. — Comparison between iV-body results and some analytical predictions in the case of 
WMAP 5 year cosmological parameters. The results at redshifts z = 1 up to k = 0.4 /iMpc -1 
are shown. We show the ratio of the predicted power spectra to the smoothed reference spectra, 
P(k)/ P nw (k), (Blue dashed, green dashed, black solid, orange solid, red solid lines, and black 
symbols are, respectively, 1-loop SPT, Regularized T-expansion, LRT, 2nd order of WH expansion 
in Eq. (|7ip and AFWH in Eq. (|68p predictions and iV-body simulation result.), and the fractional 
difference between iV-body and analytic predicted results, [-PNbody(^) — P(k)]/ P nw (k), (Blue, green, 
orange, black and red symbols are the fractional difference between iV-body and 1-loop, Regularized 
T-expansion, 2nd order of WH expansion, LRT, and AFWH. ). 

iV-body results in more quantitative ways, we also plot the fractional differences between iV-body 
simulations and the predicted power spectrum P(k) , i.e., [-PNbody(^) ~~ P{ty]/ 'P{k) (blue: iV-body 
results versus 1-loop SPT; green: regularized T-expansion; black: LRT; orange: -P^H + ^WH' rec ^ : 
AFWH in Eq. (|68])V 

The regularized T-expansion, LRT, and ^wjj + ^wh are ver y similar that we can hardly see 
any difference. Their solutions improve the overestimation of SPT, but decay at low k soon because 
of their exponential factor. On the other hand, we can find that the main difference of our result 
from the previous works with 1-loop level is the higher order of the WH expansion, P^^ir > 2). 
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Because of these terms, the AFWH in Eq. (j68|) (red solid line) does not decay and keeps the values 
well around those from iV-body simulations on BAO scales. 



6.2. Comparison with 2-loop solutions in SPT 



One merit of our interpolation is that we can directly compare our approximate solutions with 
ones of each order in the perturbation theory. In previous works, the validity of the predicted 
power spectra has been verified only by comparing with the iV-body results. However, we can 
verify the validity of our approximations, such as Eq. (|62|) and Eq. (|66p by comparing with the 
2doop solutions in SPT. The 2doop corrections are given by 



-f2ioop = P15 + P24 + -P33 + 
where each term is calculated, respectively, as 

d 3 pi d 3 p2 



f(i) 



(74) 



P 15 (z,k) 

P24(z,k) 
P33(z,k) 



30P L (z,k) 
24 



(2vr) 3 (2tt) 3 
d 3 pi d?p2 



^5(k,Pi, -pi,P2, -P2)Pl(z,Pi)Pl(z,p 2 ) 



6 



(2vr) 3 (2vr) 3 
d?p\ d 3 p2 
(2vr) 3 (2tt) 3 



F 2 (k - pi, pi)F 4 (k - pi, pi, p 2 , -p 2 )Pl(z, |k - p 1 \)P L (z,p 1 )P L (z,p 2 ), 
[F 3 (k - pi - P2,Pi,p 2 )] 2 j Pl(£, |k - pi - p 2 \)P L (z,pi)P L (z,p2). (75) 



On the other hand, we show the corresponding approximate solutions using Eq. 
as follows, 



and Eq. (jB5) 



P15 
P24 
P33 

"21oop 



[Pi 



1 



15. 



[P24 



[P33. 



a|> 2 V 2 

2_2^ 



Pi 



13) 



ap 



ap 



<k 2 a z v )P 22 , 



k 2 rr 2 

1 1 n 



2J2 



[P21oop] 



k z a. 



ap 



22- 



13 



fcV 



V P22 + 



(76) 
(77) 



Here, we have not considered the approximation of the term 
Hoop term and we can easily compute it. 



? (1) 



because it is the square of the 



In Fig. [3l we plot the correct 2-loop solutions, their approximate solutions, and their fractional 
difference, [[P] a p ~ P]/Pnw, (P = Pi5> P24, P33, and i-2k>op)> m each panel. We plot the solutions up 
to 0.2/iMpc" 1 at z = 1, because the 2doop corrections give a good result up to about these scales 
(see ([35I)). 

For Pi 5 and P24 in the top panels, the approximate solutions respectively, are overestimated 
and underestimated by about 5% at k = 0.2/iMpc _1 . On the other hand, for P33 in the bottom 
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0.1 0.2 0.1 0.2 

k[h/Mpc] k[h/M P c] 

Fig. 3. — Comparison between the approximate solutions and precise ones in the SPT 2-loop level in 
the case of WMAP 5 year cosmological parameters at z = 1. We show the ratio of the power spectra 
to the smoothed reference spectra, P/P nw , [P] a p/Pnw, and the fractional difference, [P ap — P]/P n -w 
(top left: P = P15, top right: P24, bottom left: P33, bottom right: -P2loop)- 

left panel, the approximate solution coincides very well with the precise solution within 1%. One 
may think that since there is a large difference between P15 and P24, our approximation is not 
valid. However, remember that each correction term in the perturbation theory tends to cancel 
out, resulting in small corrections. Therefore, if the approximate solution for P15 is overestimated, 
it would be natural that there is an underestimation in other solutions such as P24 to cancel out 
the overestimated solutions. As a result, for the full 2-loop corrections in the bottom right panel, 
the fractional difference becomes within 1% up to 0.2/iMpc^ 1 . 



6.3. Comparison with closure theory 



Finally, we compare with the closure theory (second Born) in (|22l ). which is one of the best 
predictions at the moment. In addition, we plot high-resolution iV-body simulations presented 
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Fig. 4. — This figure is the same as in Fig. [21 but here we compare the predictions of the AFWH 
(red solid) with those of the closure theory (purple dashed) and TV-body simulation results (green: 
high-resolution, black: low-resolution) at some redshifts (z = 0.5, 1, 2, 3). We also the fractional 
difference between the predicted power spectra and iV-body results. The red, green and purple 
symbols are respectively iV-body (low-resolution) vs. AFWH, iV-body (high-resolution) vs. AFWH, 
and iV-body (low-resolution) vs. closure theory. 



by IM). 

In Fig. [J we plot the power spectra from closure theory (purple dashed), AFWH in Eq. ([68]) 
(red solid), and iV-body results with error bars (black symbols: low resolution; green symbols: high 
resolution) at some redshifts (z = 3.0, 2.0, 1.0, 0.5). The range of plotted scales is k < 0.4hMpc -1 . 
We also plot the fractional differences (purple: iV-body result with low resolution versus closure; 
red: AFWH; green: iV-body result with high-resolution versus AFWH) 
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Fig. 5. — The left and right panels are the same as the botton left panel in Fig. 0J We compare the 
finite truncation of the WH expansion in Eq. (|8ip (red dashed) with the AFWH in Eq. (|68p (red 
line) at z = 1. In the right panel, we further plot the 5th, 6th and 7th order of the WH expansion 
up to k = l[/i/Mpc]. 

Overall, the predictions of AFWH tend to overestimate the A^-body simulations at low-A: (BAO 
scales) slightly, and then begin to underestimate at high k. The overestimation is due to the fact 
that we computed only the 1-loop level in SPT precisely. In fact, our approximate solutions are 
slightly larger than the 2-loop SPT solutions as shown Fig. [3l Therefore, to derive more precise 
prediction on BAO scales, we need to calculate up to the 2-loop level corrections. The reason for 
the underestimation is that either the expression for the high k limit would not apply perfectly 
to the range of calculation or the subleading contributions on small scales would become effective. 
We would need to compute the higher order of the WH expansion without the approximation to 
derive the precise nonlinearity in the high-A; range. Although our results certainly give slightly less 
accuracy than that of closure theory, the difference is controlled within 1% on BAO scales. 



7. Correlation function 

We compute the two-point correlation function calculated from the power spectrum in Eq. (168)) . 
which is given by 

^ . f°° k 2 dk sin(kr) , 

Usually, we are not able to compute the correlation function in SPT, because the integrand 
function diverge. For example, in the 1-loop level of SPT, the predicted solution has the scale 
dependence of k 2 PL(k) at high k, because the approximate solutions of P13 and P22 are proportional 
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Fig. 6. — Comparison between the predicted correlation functions and iV-body results (red line: 
AFWH, black dashed: linear theory, and black symbols: iV-body simulations). The results at 
some redshits (z=0.5, 1.0, 2.0, 3.0) with the range of 70 < r[Mpc//i] < 135 are shown. We 
further plot the fractional difference between the predictions of the AFWH and iV-body results, 
[£Nbod y (r)-£(r)]/£(r). 



to k 2 PL(k). Therefore, the integrand in Eq. (|78p has the scale-dependence of sin(kr) ln 2 (fc): 



k 2 sin(A;r) 



flloopO, k) oc 



k sin(kr) 



k 2 P L (k) oc sm(kr)\n 2 (k), 



(79) 



27T 2 kr ' 2tt 2 kr 

where we have used the behavior of the linear power spectrum at high-/c, Pl(^) - ► ln 2 (fe)//c 3 . 
This solution diverge at high k, and we are not able to evaluate the integration in the range of 
< k < oo. 



On the other hand, we are able to compute the correlation function for AFWH because the 
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solution has the scale dependence like the linear power spectrum at high k: 



P A F(z,k) 



8 L (z,k) + 
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oc Pl(z, k). 
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In the first line, we dropped the terms including exponential factor. The scale-dependence of <5 3 
and P22 at high-/c are proportional to k 2 5L{k) and k 2 Pi{k). As a result, Paf (red line in the 
right panel of Fig. [5]) is proportional to Pl (black dashed in Fig. [5]), and the integrand of Eq. ([78]) 
converge like the linear power spectrum. 

Furthermore, as long as we focus on the BAO scales in real space (60 < r < 140[Mpc//i]), the 
behavior of the power spectrum on small scales in Fourier space (k > 0.2— 0.4[/i/Mpc]) contributes 
very little to the result of the correlation function. Therefore, we may truncate the WH expansion 
up to the appropriate order so that the integrand function converge to zero due to the exponential 
factor and we can easily compute the correlation function: 



Paf(z, k) 



5 L (z,k) 
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(81) 



where we have truncated the WH expansion up to the sixth order. We plot the solution of Eq. (|8ip 
in Fig. [5] at z = 1, where the difference between the AFWH in Eq. (|68p (red line) and the solution 
of Eq. (|8Tj) (red dashed) up to k < 0.25[/i/Mpc] is not visible, and the solution behaves like the 
ones of closure theory at high-fc. In the right panel of Fig. [5j we further plot the solutions of the 
fifth (orange line) and seventh (orange dashed) order in the WH expansion. 

Here, we adopt the sixth order solution of the WH expansion to compute the correlation 
function in Eq. (I78p . In Fig. El we plot the analytic predictions of the correlation function, £(r) 
(red line: AFWH; black dashed: linear theory; black symbols: iV-body results), and the fractional 
difference between the predicted correlation functions from AFWH and the TV-body simulation 
results, [£Nbod y (0 - £0")]/£(r). 

Our predictions explain the displacement of the location of the BAO peaks and the smoothing 
of their amplitudes due to the non-linear effects. As a result, the fractional difference against the 
iV-body results is within 1-3%. Almost the same results are derived even for the second order 
of the WH expansion in Eq. (|7ip . because on small scales there is very little contribution to the 
correlation function around the BAO peak. This fact is also well known also in other modified 
perturbation theories (e.g., (|22l: l35l)). 
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8. Conclusion 

We have applied the WH expansion to the evolution equation of dark matter in Newtonian 
gravity. It diagrammatically corresponds to the classification of the power spectrum in which 
each order includes all of the vertex loop contributions. It is proved that the WH expansion is 
mathematically equivalent to the T-expansion approach in the Multi-Point Propagators method. 

Even if WH expansion method is physically and mathematically useful for understanding the 
non-linearity of the evolution of dark matter, the validity of the finite truncation of the expansion is 
not clear and the difficulty in the calculation will remain. To resolve these difficulties, we proposed 
a way to include the effect of all orders by assuming that the high non-linear solutions are well 
approximated by the ones in the high-A; limit. Namely we calculate only low order terms precisely 
and replace the high order solutions with the ones in the high-/c limit. 

It has been known in RPT that the matter density and velocity fluctuations of dark matter 
are exponential in the high-/c limit. We proved again this behavior in the context of SPT using the 
WH expansion by proving that the kernel functions F and G take the form of Eq. (|40p in the high-fc 
limit. Using the approximate kernel functions F and G, we proposed an appropriate interpolation 
between high-/c and low-A: solutions, and the approximate full power spectrum in Eq. (|68p . which 
approximately include the full order of SPT. 

We compared our results with some other analytic predictions (e.g., regularized T-expansion, 
LRT, SPT, and closure theory) and iV-body simulation results. Since the WH expansion is equiva- 
lent to the T-expansion and the regularized T-expansion bases on the T-expansion, we can describe 
the first order of the WH expansion, fwjj, using the regularized T-expansion in Eq. (|70p . One of 
the difference between our result and the regularized T-expansion is the manner of interpolating 
between the high-/c and \ow-k solutions, but this difference slightly affects the predicted power spec- 
trum. Another difference is that we consider the higher order of the WH expansion approximately. 
As a result, even for the 1-loop level, the predicted power spectrum in Eq. ([68]) does not decay 
due to the exponential factor as shown in Fig. [2j and results in good agreement with the iV-body 
simulation on BAO scales. 

The validity of the various modified perturbation theory (e.g., LRT, RPT, closure theory, 
. . . ) predictions is usually verified only by comparing with the iV-body simulations. However, we 
can also verify our approximation by comparing with the solutions with the SPT 21oop level. In 
Fig|3j we showed that the fractional difference between the approximate solutions and the precise 
solutions with the SPT 2-loop level is within 1 % on BAO scales (< 0.2/iMpc -1 ) for the WMAP 
five year cosmological parameters at z = 1. 

We also compared with the closure theory which is one of the best prediction at a moment, 
and the accuracy of AFWH in Eq. (|68p is comparable to or slightly less than the ones in the closure 
theory, with the fractional difference within 1% on BAO scales. 

Finally, we computed the two-point correlation function for AFWH. We can compute the 
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correlation functions because the predicted power spectrum in AFWH converges like the one from 
linear theory. Since the contributions on small scales do not affect the values of the correlation 
function, one may use Eq. (|8ip to compute the correlation function. This solution has the same 
behavior as Eq. (|68p on BAO scales and decay on small scales due to the exponential factor in Fig.[5j 
The predicted correlation functions from the AFWH agree very well with the iV-body simulation 
results, and the fractional difference is within 1 — 3%. 

We could use and apply our results to various studies of the nonlinear evolution of dark matter 
(e.g., redshift distortion effect, bias effect, and bispectrum, etc.), because our prescription is easy 
and gives good results that are comparable to closure theory, and furthermore the computational 
time is very rapid. 
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